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Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magne- 
tohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic 
turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational 
contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically 
well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the 
kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small 
and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, 
resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to 
the development of small length scales in the irrotational part of the magnetic vector potential. 
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I. INTRODUCTION 

Most astrophysical and laboratory plasmas are good con- 
ductors. This, together with high-speed flows and large length 
scales, nearly universal in the astrophysical context, makes 
for large magnetic Reynolds numbers. In the limit of infinitely 
large magnetic Reynolds number, and for domains with closed 
boundaries, total magnetic helicity is a conserved quantity. 
Here, an analogy can be drawn with mass conservation in do- 
mains whose boundaries are closed to mass flux. Furthermore, 
in open domains, the change in total mass is governed by the 
mass flux across open surfaces. In ideal magnetohydrodynam- 
ics (MHD), a similar property holds for the total magnetic 
helicity. But unlike mass, magnetic helicity depends on the 
choice of gauge. In the special case of the advective gauge, 
the magnetic helicity flux is given by the velocity times the 
magnetic helicity density 1 , making this gauge particularly in- 
teresting for studying pointwise properties of magnetic helic- 
ity. This is an important goal of this paper. 

Magnetic helicity plays an important role in many fields of 
plasma physics and astrophysics, and has applications rang- 
ing from tokamaks and other plasma confinement machines, 
to dynamo action in the Sun and the Galaxy. Our physical un- 
derstanding of the role of magnetic helicity in MHD is greatly 
aided by concepts such as Taylor relaxation 2 , selective decay 3 , 
and the inverse cascade of magnetic helicity 4 . 

Furthermore, magnetic helicity is a crucial ingredient of the 
turbulent dynamos which are believed to be the source of the 
equipartition magnetic fields in astrophysical bodies like stars 
and galaxies 5 . In all such cases the characteristic length scales 
of the dynamo generated magnetic field exceed those of the 
fluid's energy carrying scale. In dynamo theory, the forma- 
tion of such a large-scale magnetic field is typically possible 
through the a effect, which is non-zero for helical turbulent 
flows. In periodic boxes with helical turbulence, the a ef- 
fect becomes strongly quenched when the (appropriately nor- 
malized) magnetic helicity in the small-scale field (i.e., scales 
that are smaller than the energy-carrying scale of turbulent 
fluid) is comparable to the helicity in the small-scale veloc- 
ity. Conservation of magnetic helicity implies that the helicity 



in small- and large-scale fields will have comparable magni- 
tudes, so the quenching of the large-scale dynamo will oc- 
cur for weak large-scale fields. This a quenching 6 - 7 increases 
with scale separation and endures for as long as magnetic he- 
licity is nearly conserved, a resistive time that scales with the 
magnetic Reynolds number ReM = UL/rj. The quenching is 
called "catastrophic" because for the Sun ReM ~ 10 9 and the 
Galaxy ReM ~ 10 15 , and their resistive timescales are prob- 
lematically long. This rapid pre-resistive saturation of the dy- 
namo generated field poses clear difficulties in applying the- 
ory to astronomical systems, but it may be possible to allevi- 
ate the problem through magnetic helicity fluxes^. It should 
also be pointed out that problems with catastrophic quench- 
ing are often not clearly seen in present-day simulationsiS^. 
While trend lines suggest that catastrophic quenching will oc- 
cur, simulations at currently achievable, low to intermediate 
ReM and scale separation have shown significant large-scale 
fields. 

There exists reasonable observational evidence is support of 
such fluxes of magnetic helicity. The Sun's surface magnetic 
field shows helical structure s 13,14 . Further, it was showni^ that 
the S-shaped (helical) regions which are active in the corona 
are precursors of coronal mass ejections (CMEs) and later 16 
that those regions are more likely to erupt. This suggests that 
the Sun sheds magnetic helicity via CMEs. Since the Sun's 
large-scale magnetic field is believed to be generated by a 
helical dynamo 17,18 this shedding of magnetic helicity could 
play an important role in the 1 1 year solar cycle. Physically, 
magnetic helicity fluxes out of the domain can be mediated in 
many ways, such as the aforementioned CMEs for the Sun 19 
or fountain flows in the case of galaxies^ 7 -. In direct sim- 
ulations magnetic helicity fluxes are permitted by adjusting 
the boundary conditions, e.g., to vertical field boundaries, but 
their actual presence can be difficult to ascertain. Internal he- 
licity fluxes have also been found to alleviate a quenching 18 
in systems with internal boundaries that separate zones of op- 
positely signed kinetic and magnetic helicities. 

A difficulty in addressing the generation and transport of 
magnetic helicity is its gauge dependence. We denote the 
magnetic vector potential as A such that B = V x A is 
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the magnetic field. Magnetic helicity H = J v A ■ B dV is At the core of MHD is the induction equation, 
independent of the gauge for perfectly conducting boundaries, 
as well as periodic boundaries so long as A is also required 
to be periodic. However, if one wishes to study the transport 
of magnetic helicity for physically motived systems a non- 
volume integral formulation will be needed. Magnetic helic- 
ity density, h = A ■ B, the quantity we will be working with, 
clearly depends on the gauge choice for A. The gauge de- 
pendence of fluxes of mean magnetic helicity contained in the 
fluctuating fields was examined via direct numerical simula- 
tions (DNS) for three different gauges 20 , and it was found that, 
averaged over time, they do not depend on the gauge choice. 
This is a result of the fact that, for sufficient scale separation, 
the magnetic helicity of the fluctuating field can be expressed 
as the density of linkages, which in turn is gauge-invariant 21 . 
This result implies that the study of specific but useful gauge 
choices is a meaningful task. 



(1) 



where U is the velocity and rj is the molecular magnetic dif- 
fusivity. Equation ([T]) can be uncurled to give an evolution 
equation for the magnetic vector potential A, but only up to 
a gauge choice. In the Weyl gauge, indicated by a superscript 
W on the magnetic vector potential, we just have 



dA w 

dt 



= U x B i]J, 



(2) 



but by adding the gradient of a scalar field, the vector potential 
can be obtained in any other gauge. Of particular interest to 
this paper is the advective gauge, 



In this work we examine the properties of magnetic helic- 
ity density in a particularly interesting gauge-family which 
we call "advective" because in this gauge the effect of ve- 
locity on the evolution equation of magnetic helicity takes 
the form of a purely advective term. In previous work-i this 
gauge choice was shown to be crucial to understanding mag- 
netic helicity fluxes in the presence of shear, including the 
Vishniac-Cho flux 22 . Unfortunately, evolving A in this gauge 
proves numerically unstable. This may be related to earlier 
findings in smoothed particle MHD calculations 23,24 . There, 
the problem was identified as the result of an unconstrained 
evolution of vector potential components, which was argued 
to be connected with "poor accuracy with respect to 'reverse- 
advection'-type terms" 23 . Our present work clarifies that this 
instability is related to the excessive build-up of irrotational 
contributions to the magnetic vector potential. These contri- 
butions have no physical meaning, but discretization errors at 
small length scales can spoil the solution dramatically. 

We shall therefore describe a novel method for obtaining A 
in this gauge by evolving it first in a numerically robust gauge 
and then applying a gauge transformation with a simultane- 
ously evolved gauge potential. This will be referred to as the 
A method throughout the text. Next, we show that the mag- 
netic helicity density in the advective gauge tends to be small 
even pointwise, provided turbulent effects are still weak, and 
discuss the analogy with passive scalar transport. We con- 
clude by pointing out that resistive terms break the analogy 
with passive scalar advection through the emergence of a tur- 
bulently diffusive magnetic helicity flux. 



II. MAGNETIC EVOLUTION EQUATIONS 

A. Weyl and advective gauges 

In this work we remain within non-relativistic MHD and 
hence neglect the Faraday displacement current. So the cur- 
rent density is given by J = V x B, where B is the magnetic 
field and we use units where the vacuum permeability is unity. 
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where A W:a is the gauge potential that transforms from A w 
to A a . We demand that 25 
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Here, B/Bt — d/dt + U ■ V is the advective derivative. Con- 
sequently one can show that A W:a obeys the evolution equa- 
tion (see Appendix lAl 
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Thus, to obtain A a , one can either solve Eq. © directly or, al- 
ternatively, solve Eq. (f2j) together with Eq. © and use Eq. (0 
to obtain A a . A possible initial condition for A W:a would be 
A W:a = 0, in which case A a — A w initially. For numer- 
ical reasons that will be discussed in more detail below, we 
shall consider the indirect method of obtaining the magnetic 
vector potential in the advective gauge, but starting from more 
numerically stable gauge which will be discussed in § 111 Bl 

Variants on the advective gauge have seen significant use, 
particularly in DNS with constant imposed shear. Although 
the magnetic field in such simulations must obey shearing- 
periodic boundary condition the vector potential need not. 
In particular, the evolution equation (J2j does not impose 
shearing-periodicity on the vector potential, while Eq. (0]) 
does, enabling shearing-periodic numerical simulations^ in 
terms of A. 

For our purposes, the importance of Eq. (|4]i lies in the form 
of the magnetic helicity density evolution equation. By writ- 
ing the induction equation in the form 
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computing B{A a -B)/Bt = A a -BB/Bt + B- BA a /Bt, and 
noting that the AiUi.jBj terms from both equations cancel, 
we find that 



Bh a 
~Bt 



-/i a V • U - V ■ (r)J x A a ) - 2r/J ■ B, 



(7) 
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which shows that in ideal MHD (77 = 0) under the assumption 
of incompressibility ( V • U = 0) the magnetic helicity density 
in the advective gauge, h a — A a ■ B is just advected with the 
flow like a passive scalar, i.e. 



Dh a 



(for 77 = and V- (7 = 0). 



(8) 



In the general case with V • U ^ 0, the rate of change of the 
local value of h a is given by — V • (h a U), which is analogous 
to the continuity equation for the fluid density. However, for 
1) 7^ 0, there is also a source term, 



dh a 

^- = -2r,J-B-V-F a 

at 



(9) 



as well as a resistive contribution to the magnetic helicity flux, 



F a = h a U + r]J x A a . 



(10) 



In this paper we address the question how the 77 J x A a contri- 
bution scales in the limit 77 — > 0, i.e. for large values of ReM- 
It could either stay finite, just like the resistive energy dissipa- 
tion 77 J 2 , which tends to a finite limit 5 as 77 — »• 0, or it could 
go to zero like the source term r]J ■ i ? 27,28 . 



and finally using the gauge transformation Eq. ( TT3b . For a full 
derivation of this equation we refer to Appendix iBl Note that 
the microscopic magnetic diffusivity automatically enters the 
A r:ar equation as a diffusion term, which implies that the A r:ar 
equation is numerically well behaved. 

The magnetic helicity density h al = A al • B in the advecto- 
resistive gauge can be calculated from the magnetic helicity in 
the resistive gauge through h al = /i r +VA r:ar -B, and it obeys 



dh a 
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-2r)J B - V -F a 



with 



h ai 'U - r)(V ■ A ai )B + r)J x A a 



(16) 
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For comparison, the evolution equation of the magnetic he- 
licity density in the resistive gauge is given by an equation 
similar to ( TToT ). but with h al being replaced by h T and F ar 
being replaced by 

F 1 = h r U - (U ■ A 1 ' + 77V ■ A*)B + 77J x A 1 , (18) 

which contains a non-advective velocity driven flux of the 
form (U ■ A r )B - even in the ideal case. 



B. Resistive and advecto-resistive gauges 

There are two important issues to be noted about the equa- 
tions discussed above. Firstly, for numerical reasons, Eq. (f2]i 
is often replaced by 



dA r 
~dt 



= U x B + TiV 2 A r , 
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where A 1 ' is the magnetic vector potential in the resistive 
gauge and we have assumed that 77 = const; otherwise there 
would be an additional gradient term of the magnetic diffusiv- 
ity that results from22 



rjJ + V(i]V ■ A) = r/V 2 A + (V • A)V?7. 



(12) 



This "resistive" gauge introduces an explicit, numerically sta- 
bilizing diffusion term for each component of A. Secondly, 
and again for numerical reasons, Eq. (O should be solved with 
a small diffusion term proportional to V 2 A W:a . These two is- 
sues are actually connected and can be resolved by consider- 
ing the gauge transformation 



A 1 + VA r 



(13) 



which allows us to obtain the magnetic vector potential A a 
in the advecto-resistive gauge obeying 
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-^ = - Uj , i Af + r ) W 2 Ar, 
by solving Eq. ( fTTT i for A 1 together with 
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C. Numerical details 

We perform simulations for isotropically forced, triply pe- 
riodic cubic domains with sides of length 2ir, as was done in 
earlier work 28 . The 77 J • B term in ((9} implies (and past simu- 
lations have shown) that such a system will experience a slow, 
but steady production of magnetic helicity. This is the price 
to pay for a system which is both helical, providing us with 
a signal, and homogeneous, so avoiding extraneous magnetic 
helicity fluxes. In addition to the uncurled induction equation 
( fTTT ) and the gauge transformation evolution equation ( fT5b . we 
solve 
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where c s (= const) is the isothermal sound speed, p is the 
density, .F v i sc = p _1 V • {2pvS) is the viscous force, S,j = 
+ Ujj) — ^Sij'V ■ U is the rate of strain tensor, v is the 
kinematic viscosity, / the forcing term, and cl = 1 is a pref- 
actor that can be put to to turn off the Lorentz force in kine- 
matic calculations. As in earlier work 28 the forcing function 
consists of plane polarized waves whose direction and phase 
change randomly from one time step to the next. The modu- 
lus of its wavevectors is taken from a band of wavenumbers 
around a given average wavenumber fcf. The magnetic vec- 
tor potential is initialized with a weak non-helical sine wave 
along one direction. In some cases we shall also consider so- 
lutions to the passive scalar equation in the incompressible 
case, 
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where k is the passive scalar diffusivity. Following earlier 
work2£, we impose a linear gradient in C, i.e. C = Gz + c, 
and solve for the departure from this gradient G, i.e. 



— = kV 2 c-GU z , 
Dt 



(22) 



where GU Z acts essentially as a forcing term. 

We use the Pencil Code 

(http://pencil-code.googlecode.com: 31 to solve the equa- 
tions for A 1 , U, A r:ar , p, and in some cases also c. The 
calculations involving A r:ar have been carried out with 
the publicly available revision |rl521l! (or similar) of the 
module |special/ advective_gauge . f 90| 

The control parameters we use are the magnetic Reynolds 
number ReM, the magnetic Prandtl number PrM, and the 
Schmidt number, 



ReM = — j—, Pr M = -, Scee-, 
•qkf r\ k 



(23) 



where u rms is the root mean square velocity. We use fcf = 3/ci 
where fci, the box wavenumber, is unity. The numerical res- 
olution is varied between 32 3 and 256 3 meshpoints for val- 
ues of Re and ReM between 3 and 300. In one case we used 
ReM ~ 800, which was only possible because in that case 
we used PrM = 10, so that most of the energy gets dissi- 
pated viscously, leaving relatively little magnetic energy at 
high wavenumbers 3 ^. 



III. IMPORTANCE OF MAGNETIC HELICITY DENSITY 

A. Implications of for dynamo theory 

Magnetic helicity is not only of interest by being a con- 
served quantity in ideal MHD, but also by being the basis 
of a methodology to treat nonlinear helical MHD dynamos, 
namely dynamical a quenching 33 . This methodology relates 
the current helicity in small scale fields with the magnetic he- 
licity in small-scale fields, j ■ b ~ k 2 a ■ b, and invokes the 
magnetic a effect 4 . The evolution equation of the magnetic 
helicity density then becomes the evolution equation of the 
magnetic part of the a effect and the nonlinear evolution of 
the dynamo can be modeled. This methodology has been used 
successfully in systems where no net helicity flux is possible, 
and initial work invoking the methodology has captured the 
behavior of at least one system with finite helicity fluxes 3 ^. A 
major prediction of the theory is that in the absence of pref- 
erential helicity fluxes of small-scale fields, dynamo action is 
quenched to sub-equipartition mean field strengths. This phe- 
nomenon is sometimes referred to as "catastrophic quench- 
ing". 



fluxes, Eqs. ( TTOb and ( TPTI i . In the limit of ideal, incompress- 
ible, kinematic MHD, Eq. (0 is the evolution equation for a 
passive scalar. Even in non-ideal MHD, if the fluctuations of 
/i ar due to the velocity field U were purely advective in nature 
(i.e. passive), magnetic helicity transport would only be resis- 
tive, large-scale advective, and/or turbulently diffusive. This 
would forbid the preferential export of small-scale magnetic 
helicity and might call for alternate solutions to the catas- 
trophic quenching problem than helicity fluxes^. 

While in ideal MHD (rj — 0) the resistive terms in (0 van- 
ish, resistive terms need not vanish in the limit of rj — > 
(high ReM). For example, in a turbulent flow, Ohmic dissi- 
pation rjj 2 tends to a finite value as -q decreases. The need 
for non-resistive solutions to the build-up of magnetic helicity 
is therefore not a given. We will examine this by performing 
kinematic simulations where the Lorentz force is turned off, 
i.e. cl = 0. 

If the Lorentz force is significant, the fluctuations of /i al and 
U might be correlated beyond simple turbulent diffusion con- 
cerns (i.e. the fluctuations of h m ' could drive flow patterns). 
In the limit of incompressible flows, if the helicity is uniform, 
then the only source terms for helicity patterns of finite k are 
the resistive terms. The terms are small compared to dimen- 
sional estimates for the velocity terms when ReM S> 1. We 
will look for signals of magnetic helicity transport by exam- 
ining spectra of h r and /i ar as (pseudo) scalars, together with 
spectra of a true passive scalar. As we will show, the advecto- 
resistive gauge is adequately efficient at turbulently diffusing 
magnetic helicity that no inertial range for the magnetic he- 
licity density can be identified. However, the spectra of h T 
help elucidate previous results 34 which found diffusive fluxes, 
but at values well below turbulent diffusivities. Instead, our 
spectra show clear diffusive behavior in the inertial range, but 
the mere existence of the inertial range implies non-diffusive 
behavior. 

We emphasize that our spectra of h r and h m ' have nothing 
to do with the usual magnetic helicity spectrum that obeys a 
realizability condition and whose integral gives the volume- 
averaged magnetic helicity. Here we are looking instead at 
the power of the magnetic helicity density as a (pseudo) scalar 
field. Our hk measures the spatial variation of h. In order 
to avoid confusion, we shall refer to these spectra as scalar 
spectra. 



IV. RESULTS 

The results reported below for the magnetic helicity density 
h refer to the advecto-resistive gauge and have been obtained 
by the A method, unless indicated otherwise. The results from 
the direct method agree (£> |IV At , but this method develops an 
instability when nonlinear effects become important (§ |IV B| >. 



B. Magnetic helicity as passive scalar 

In the advective and advecto-resistive gauges, the veloc- 
ity appears in the evolution equations of the magnetic helic- 
ity density, Eqs. (O and ( TTol l. only as advection terms in the 



A. Agreement between A and direct methods 

To test the agreement between the A method and directly 
solving the induction equation in the advecto-resistive gauge, 
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FIG. 1: (Color online) Time dependence of the normalized helic- 
ity for the advecto-resistive gauge with the direct method and the A 
method. Both curves agree perfectly just until the moment when the 
code develops an instability in the direct calculation. Time is normal- 
ized in terms of the magnetic diffusion time. The fit is an exponential 
relaxation to a constant value proportional to 1 — exp(—2r]k% n lS.t), 
where At = t — t S at is the time after the small-scale magnetic field 
has saturated^ and fc m = 1.4fci has been chosen for a good fit. 




we plot the normalized rms magnetic helicity h^ s with re- 
spect to time (Fig. [TJ. Note that the non-dimensional ratio 
fcl^rms/^rms nas a well-defined plateau during the kinematic 
stage. Below we shall study the average value of this plateau 
as a function of magnetic Reynolds and Prandtl numbers. At 
the end of the kinematic phase, there is a slow saturation phase 
on a resistive time scale during which the large-scale field of 
the dynamo develops 28 . The results of the two calculations 
agree just until the moment when the direct calculation devel- 
ops a numerical instability, whose nature will be discussed in 
more detail below. The perfect agreement until this moment 
can be taken as confirmation that the A method works and is 
correctly implemented in the code. 



B. Nature of the instability 

In Fig.|2]we show time series for a range of modest values 
of ReM and two resolutions, 32 3 and 64 3 . Reducing the mag- 
netic Reynolds number may stabilize the system somewhat, 
but changing the resolution has no clear effect. In Fig. [3] we 
present data from equivalent runs that solve either (fl4] > or al- 
ternatively ( ITTb and (15[ . We can see that the solutions match 
up until time t = 220/c s fcf, where the run that solves ( fl~4l 
becomes unstable. 

The key point is that when we evolve dTTb and ( fl~5l l, A never 
enters the equations for physical quantities. However, when 
we evolve (fl~4l >. the magnetic field includes a term V x (VA), 
which, when computed numerically, is not zero. The first 
panel in Fig. [3] shows the power spectra of the vector poten- 
tial. Comparing the advecto-resistive gauge (dashed/red) with 



FIG. 2: (Color online) Evolution of _B rms /_B cq for small values of 
ReM between 4.3 (top) and 2.1 (bottom), using 32 3 (solid lines) and 
64 3 (dashed, red lines). In each case, time on the abscissa is normal- 
ized by the growth rate A, whose value is given in each panel in units 
of the inverse turnover time, r~ = u lnlB kf. The ends of each line 
mark the point when the solution became unstable. 

resistive gauge (dotted/blue) we see that A ar = A r + VA has 
significantly more power at high k than A 1 . Numerics can- 
not adequately handle the requirement that V x VA = at 
high k in the direct method, introducing errors in B, as can be 
seen in the second panel. This fictitious increase in magnetic 
power at high k (and the attendant increase in current) result 
in a fictitious high k increase in the velocity field (third panel) 
that produces the numerical instability. The results of Fig. [2] 
suggest that the power of A (remembering that J includes that 
the third derivative of A) drops slowly enough at high k that 
numerical stability can only be achieved by enforcing an ade- 
quate resistivity r\ to damp A for only modest wavenumbers. 
Indeed, any gauge with large power in A for high k is ex- 
pected to be numerically unstable, and the method sketched 
in Appendix lAl or IB"! may be used to make the connection be- 
tween analytical results in such a numerically unstable gauge 
and numerical results produced in a stable gauge. 

C. Evolution of rms helicity density 

In Fig. [4] we present a time series of the normalized rms 
magnetic helicity density in the kinematic regime (Lorentz 
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FIG. 3: (Color online) Power spectra of A, B, and U for two runs 
that are identical except that the first run solves for A al directly while 
the second solves for A 1 and A. In the top panel we plot the spec- 
trum of A obtained either via A al — A 1 + VA (dashed) or di- 
rectly, A^ ir (solid/red), and compare with A 1 (thick gray/yellow), 
showing that the vector potential in the advecto-resistive gauge has 
much more power at high k. The inset shows the time evolution of 
the normalized h m s shortly before the time of the numerical insta- 
bility. The dash-dotted line indicates the time for which the power 
spectra are taken. In the second panel we present magnetic energy 
spectra obtained in the direct gauge (solid/red), with the A method 
(dashed/black) as well as fcA ar (dotted/blue), showing that there is 
significant power in the irrotational part of A. We see that in the 
direct calculation of A 3,1 ' the numerics is unable to adequately han- 
dle the high wavenumber power of A al with consequences for the 
velocity seen in the last panel (solid/red line). The spectra of B and 
U agree for resistive and advecto-resistive gauges (thick gray/yellow 
line underneath the dashed black line) because the evaluation of the 
curl of a gradient has been avoided (last two panels). The three spec- 
tra are all taken for t = 210/c a fei. 



force turned off, i.e. cl = 0). In both the advecto-resistive 
and resistive gauges, there is an initial adjustment of the non- 
dimensional ratio k\ h rms / B^ ms to a certain value, followed by 
a plateau. In the kinematic regime the magnetic helicity den- 
sity is passive and the advection term in the advecto-resistive 
gauge merely serves to turbulently diffuse any local concen- 
trations of h ar . Therefore there cannot be any spontaneous 
growth of h aT , except for effects from the resistive terms in 
the early adjustment phase. Turbulent diffusion itself, on the 
other hand, cannot generate variance of h ar . 



FIG. 4: (Color online) Time dependence of the rms values for the he- 
licity in the advecto-resistive (solid/red) and resistive (dashed/blue) 
gauges with the Lorenz force switched off, i.e. cl = in both cases. 



In Figs.|5]and|6]we plot the height of the rms-magnetic he- 
licity density plateau as a function of ReM for several values 
of the magnetic Prandtl number and constant forcing ampli- 
tude. The differences between the evolution equations for h r 
and h al are contained entirely in the flux terms so the volume 
integral of h is the same in the two gauges. Any difference 
between the rms values of h therefore is due to spatial fluctu- 
ations generated by the flux terms. 

We fit the data points in Fig. [5] with functions of the form 
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The fit results for the parameters are presented in Table U Of 
interest is c, which increases with PrM and scales approxi- 

1 12 

mately with Pr^ . A more general, albeit less accurate fit is 
given by 



fci/i a 



3 Re 



M 



Re M /PrM 3 
50 



(25) 



see Fig. [7] 

It is clear that high wavenumber fluid eddies (which are 
damped for small Re, i.e. large Pi"m, contribute significantly 
to hf ms for ReM > 100, while from Fig. [6] we see that they 
do not contribute to ftf ms . That these eddies could contribute 
in the advecto-resistive gauge is to be expected as the advec- 
tive nature of that gauge implies the existence of an efficient 
turbulent cascade; the fact that they do contribute there and 
that the r/J x A 3,1 ' and t/(V • A ar )B terms remain important 
implies that resistive terms both become important at small 
length scales and have non-dissipative effects. This is ex- 
plained by the fact that A m ' develops a strong high-fc tail; see 
also Fig. [3] This is confirmed in Fig. [8] which shows that 
the resistive magnetic helicity fluxes in the advecto-resistive 
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FIG. 5: (Color online) ReM dependence of fci/irms/ -Brms f° r the 
kinematic phase. Values are averages over times where they reach a 
stationary state. The curves represent fits according to l l24l >. 



FIG. 6: (Color online) ReM dependence of k\h T TV0S /B'} ms for the 
kinematic phase. Values are averages over times where they reach a 
stationary state. A —1/4 power law can be seen. 



gauge are proportional to ReM- In this gauge the rms resistive 
helicity fluxes are therefore independent of the actual value of 
the resistivity, staying finite even in the high ReM limit. This 
is quite different from the resistive magnetic helicity fluxes in 
the resistive gauge, and the global magnetic helicity dissipa- 
tion (which is gauge-independent): both terms are only pro- 
1 1% 

portional to Re^ and, after multiplying with r\ these terms 
tend to zero for ReM — > oo. 



D. Comparison with passive scalar 

In Fig. [9] we present scalar spectra of the magnetic helic- 
ity density for both the resistive and advecto-resistive gauges 
and for the passive scalar concentration c, in the kinematic 
(arbitrary units) and saturated regimes. The passive scalar 
spectrum shows a peak at the forcing scale, k{/k% = 3, fol- 
lowed by an approximate fc -5 / 3 subrange and an exponential 
diffusive subrange. As long as the magnetic energy density 
is still small compared with the kinetic energy density, the 
field exhibits exponential growth and a Kazantsev fc 3 / 2 energy 
spectrum, which is well seen in simulations even at magnetic 
Prandtl numbers of unity both with and without kinetic helic- 
ity in the velocity field 35 . This fc 3 / 2 spectrum is also reflected 
in the scalar spectrum of /i ar . The scalar spectrum of h T is 



TABLE I: Fit parameters for equation i24\ and Fig. [5] 



Pr M 


a 


b 


c 


line type 


1 


0.7 


3 x 10" 3 


1.2 


solid/blue 


5 


0.9 


4 x 10~ 4 


2.0 


dashed/green 


10 


1.0 


5 x IO -5 


3.5 


dotted/red 




1 10 100 

Re M /Pr> /3 



FIG. 7: (Color online) Dependence of fci farms/ B^ ms , scaled by Pr^ 3 

1/3 

on ReM /Pi"m f° r me kinematic phase and PrM = 1 (filled circles), 
5 (open circles), and 10 (plus signs). The solid line represents the fit 
ofEq. {25). 



somewhat steeper and closer to k 2 , indicating that h 1 is dom- 
inated by white noise in space at large scales. 

The saturated regime exhibits some interesting properties. 
The pronounced peak of the power of the passive scalar at the 
driving scale is easily understood as being due to the source 
of c. However, the magnetic helicity density in the resistive 
gauge shows a significant peak there as well, while it does not 
in the advecto-resistive gauge. This implies that the veloc- 
ity term in Eq. ( fl~8l ) generates significant spatial variations in 
the magnetic helicity density - even in the absence of external 
modulations. As in dynamical a quenching, h influences the 
a effect, this suggests a way to quantify the appropriateness 
of different gauge choices: systems where spatial and tempo- 
ral fluctuations in a can be adequately constrained would al- 
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FIG. 8: (Color online) ReM scaling of the rms value of J x A, nor- 
malized by ReM5 r 2 ms , for the advecto-resistive and resistive gauges. 
The solid line represents constant scaling, i.e. r\J x A aT ~ const, 
while the dashed line represents inverse square root scaling, i.e. 

— 1/2 

r\J x A 1 ' oc Re M , for three runs with PrM = 1 in the saturated 
regime. The dotted/blue line shows that rjj 2 , properly normalized, 
is approximately constant. 
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low one to determine whether spatial fluctuations in h, as seen 
in Fig. [9] are fictitious as suggested by the advecto-resistive 
gauge or not. 

The spectra of /i ar in the saturated regime does not present 
a clear inertial range, so we cannot draw strong conclusions 
as to possible non-diffusive turbulent fluxes. However, h r fol- 
lows the same cascade as the passive scalar. Previous studies 
in that gauge 1 found that magnetic helicity fluxes were best 
treated as diffusive, although the fits were imperfect. The dif- 
fusive nature is clearly seen in the spectrum while the imper- 
fections of the diffusive fit can be seen in the generation of a 
peak at the driving scale. This evidence in support of diffusive 
magnetic helicity fluxes gives us the confidence to predict at 
what ReM diffusive magnetic helicity fluxes will play a dom- 
inant role in dynamo saturation, i.e. when the diffusive fluxes 
have a greater effect on magnetic helicity evolution than the 
resistive terms. This will be done in §[V]where we re-analyze 
simulation data from earlier work^i. 



V. REVISITING EARLIER WORK 

Earlier wor k 20 ' 34 on magnetic helicity fluxes in inhomoge- 
neous open systems confirmed that the magnetic helicity den- 
sity of the small-scale field is gauge-invariant - even if that of 
the large-scale field is not. The divergence of the mean mag- 
netic helicity flux of the small-scale field is then also gauge- 
invariant, but its value is small compared with resistive mag- 
netic helicity dissipation. We return to this work to estimate 
at what ReM diffusive magnetic helicity fluxes will begin to 
play a dominant role in dynamo saturation. 

We emphasize that we are now discussing helicity prop- 
erties of what we call the small-scale field. Such a field is 



FIG. 9: (Color online) Power spectra of h 1 , /i ar , and the passive 
scalar c, both in the kinematic regime (top) and the nonlinear satu- 
rated regime (bottom) for Re = 80 with PrM = Sc = 1. In the kine- 
matic regime, the dash-dotted lines have slopes +2 for h T , +3/2 for 
h ar , and —3/2 for c (top) and —5/3 for c in the saturated regime. 

defined by introducing an averaged magnetic field, B, indi- 
cated by an overbar. Following earlier work 20,34 we restrict 
ourselves here to planar (or horizontal) averaging. The small- 
scale field is then given by b — B — B, and the mean mag- 
netic and current helicity densities of the fluctuating fields are 
then hf = a b and j ■ b, respectively, where V x a = b 
and j = V x b. Turbulent diffusion and the a effect imply 
helicity transfer between scales 36,37 through the mean electro- 
motive force of the fluctuating field, £ = u x b, so that the 
evolution equation for hf takes the form 

—f = -2£-B-2T 1 Tb-V-F { . (26) 
at 

Here, both hf and V • Ff are a gauge-dependent, but if there 
is a steady state, and if h{ is constant, then dhf/dt = 0, and 
since both £ • B and j ■ b are gauge-invariant, V • Ff must 
also be gauge-invariant. Numerical values for £ ■ B, j ■ b, and 
V • Ff were given earlier 34 for a particular simulation of a slab 
of helically driven turbulence embedded in a poorly conduct- 
ing non-helically driven turbulent halo. In Fig.[l0]we show 
the scaling of all three terms versus ReM- Note that —£ ■ B 
is balanced mainly by j ■ b. However, if the current trend, 

i — 1/2 

j ■ b ~ Re M and V • Ff ~ Re M were to continue, one 
might expect a cross-over at ReM ~ 3 x 10 4 . If so, the scal- 
ing of £ ■ B is expected to become shallower, following that 
of V ■ Ff . Given that the largest ReM accessible today is of 
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FIG. 10: (Color online) Scaling of £ B, j ■ b, and V-ff versus ReM 
for the data of an earlier simulation 34 of helically driven turbulence 
embedded in a poorly conducting non-helically driven turbulent halo. 
The symbols show actual data obtained from simulations, the dashed 
lines are the extrapolation to high ReM- 

order 10 3 , we may conclude that an alleviation of quenching 
through diffusive magnetic helicity fluxes will not be promi- 
nent in simulations for the near future. Nevertheless, astro- 
physical systems such as the Sun are orders of magnitude be- 
yond the estimated critical point of ReM ~ 3 x 10 4 ; and we 
expect their dynamo dynamics to behave accordingly. 

VI. CONCLUSIONS 



for the gauge transformation. 

The advecto-resistive gauge has allowed us to examine both 
the consequences of finite resistivity for magnetic helicity 
density as well as the possibilities of turbulent transport. The 
magnetic helicity flux, and in particular the contribution from 
7] J x A al (properly normalized) reaches a constant value as 
r] —> 0. This behavior is similar to the behavior of energy 
dissipation in turbulence, known as the law of finite energy 
dissipation 38 . This is interesting as the source term for the 
volume integrated magnetic helicity H does in fact tend to 
zero as T) does. In this sense, the high ReM behavior of mag- 
netic helicity is richer than previously anticipated. Indeed, 
the generation of spatial magnetic helicity fluctuations ex ni- 
hilo in non-advecto-resistive gauges is interesting, with po- 
tentially testable implications. We expect that the magnetic 
helicity fluxes resulting from terms of the form 77 J x A ar can 
be modeled as turbulent Fickian diffusion-type fluxes down 
the gradient of mean magnetic helicity. However, it is clear 
that fluxes from turbulent diffusion provide only a poor es- 
cape from catastrophic a quenching, partly because they can- 
not distinguish between large- and small-scale fields. Further- 
more, in simulations with such turbulent diffusion fluxes, their 
contribution is still much smaller than the local resistive mag- 
netic helicity dissipation 2034 . However, the latter decreases 
faster (~ Re^ 1 ) with magnetic Reynolds number than the for- 

— 1/2 

mer (~ Re M ), so one may estimate that only for magnetic 
Reynolds numbers of around 10 4 one has a chance to see the 
effects of turbulent diffusion. If true, however, such fluxes 
would definitely be important for the magnetic Reynolds num- 
bers relevant to stars and galaxies - even though such values 
cannot be reached with present day computer power. 



In view of the fact that the time averaged magnetic helicity 
of the fluctuating fields is gauge-invariant in systems with suf- 
ficient scale separation, the gauge-freedom can be exploited to 
gain insights using gauges that are particularly revealing. Here 
we have examined an interesting gauge, the advecto-resistive 
gauge. As the advecto-resistive gauge is inherently numeri- 
cally unstable, we had to implement a possibly universal tech- 
nique to run numerical simulations in such unstable gauges by 
running in a stable gauge while also solving a further equation 
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Appendix A: Derivation of Eq. 

We begin by expressing U x B in terms of A, 

(UxB) i = U j A j , i -U j A i j. (Al) 

The last term can be subsumed into an advective derivative 
term for A. Using furthermore UjAj^ = (UjAj).i — Uj^Aj, 
we can write Eq. (O as 

DA W 

—^^-U J , l AY + (U-A w )^-r h J l . (A2) 
We now insert Eq. ® for A w = A a - VA W:a , so 



. W:a 



DA a DAy :a 

Dt Dt J ' 3 h > 3 

+{U-A yr ) ti -r i J i . (A3) 



and note that 

DA W:a 



Dt 



|^ DA W:a ^ 



(A4) 



The last term cancels and we are left with 



DA 



< DA W:a 



D/ , U jti Aj + r)Ji = V, ^-^f- + U ■ A w j , (A5) 

so we recover the evolution equation for the advective gauge 
provided Eq. (0 is obeyed. 



Appendix B: Derivation of Eq. dTH 

We present here the derivation of the transformation from 
the resistive gauge to the advecto-resistive gauge, proceeding 
analogously to the derivation presented in AppendixlAl How- 
ever, instead of Eq. (IA2I) we now have 



Dt 



U ]il A) + (U-A"), l + r 1 V 2 ^ 



(Bl) 
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Inserting Eq. ( H~3T > for A 1 = A ar - VA r:ar , we obtain an 
Equation similar to jA3i , 



DA? 
Dt 



DA T l a 
Dt 



-U jti Af + r.,.V;' : +(U-A I 



+r 1 V 2 A? - rjV z A 



2 a r:ar 



(B2) 



so we recover the evolution equation for the advecto-resistive 
gauge provided Eq. (15[ is obeyed. 



which leads to 

DA? 1 



2 A&T 



Dt 
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2 a r:ar 
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(B3) 



